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Analytical expressions are provided for the configurations of an inextensible, flexible, twistable 
inertial string rotating rigidly about a flxed axis. Solutions with trivial radial dependence are helices 
of arbitrary radius and pitch. Non-helical solutions are governed by a cubic equation whose roots 
delimit permissible values of the squared radial coordinate. Only curves coplanar with the axis of 
rotation make contact with it. 



The classical problems of thin strings have invited attempts at solution since at least the seventeenth century. 
They include the catenary, velaria, lintearia, and motion in a radial potential [1 . Closely related to the last are 
the rigid rotations of a string about a fixed axis. This problem and its generalizations, including additions of a 
uniform gravitational force, air drag, end masses, and tangential motion of material elements along the string, have 
been considered in diverse contexts by engineers, physicists, and mathematicians. There have been abstract studies 
^irS , as well as models of textile manufacture [9lll7j. payload manipulation by space-, air-, and marine craft |l8n25] . 
jump (skipping) ropes and wind turbines [^BHSO] . Whether thought of as tethers whirling drogues, yarn balloons, or 
troposkeins, the shapes are generated by a simple and ubiquitous process. Thus, it is rather surprising that a fully 
three-dimensional analytical exploration of even the bare-bones problem remains unpublished. This situation is likely 
due to the imposition of certain boundary conditions that severely restricted the solution space explored by previous 
researchers. This note details solutions of the unaugmented problem. 

Given a time-dependent curve X(s,t) parametrized by arc length s, the balance of inertia and line tension in an 
inextensible, perfectly flexible and twistable string of uniform mass density /i is described by the vector wave equation 
and metrical constraint 
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where the stress cr is a multiplier field enforcing the constraint [71 l31f[57] . I will use two frames to describe the curve. 
The first is the moving triad of unit vectors (t, n, b), defined such that 
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The two extrinsic curvatures k and t will not be important in what follows. Velocities in this frame are represented 
by ax = Tt + iVn + Bh, with inextensibility implying dsT = kN. The second frame is the Cartesian (x, y,z), in 
which I will represent the curve as 



r{s) cos[wt + (j){s)] 
X = I r{s) sm[ujt + 0(s)] 



(4) 



a shape consisting of material elements rotating around the z axis with nonzero angular velocity uj, with no overall 
fiux of material along the string. In this representation, the constraint equation ([2| is (a?")^ + (a^)^ + {rds(f>)^ — 1. 
Finally, I will denote any constant of integration by a subscripted c. 



with constant R and P, are solutions 



Before proceeding, note that all helices, r — R, (f> = Ps, z = ± (l 

of ([T]) and ([2]) bearing a uniform stress cr = independent of radius. This includes circles perpendicular to, and 
straight lines parallel to, the axis of rotation. The latter have infinite stress, a pathological limit akin to the straight 
catenary. 

Now consider the solutions with nontrivial radial derivatives dgr. Projection of ([T]) along z immediately gives 
adgZ = ci. So if the axial slope is zero anywhere, it is zero everywhere. Noting that dtT = dfN = dtB = 0, projection 
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of ([T]) along b gives Tdti ■ b = —Ndtn ■ b. Inserting into the projections of ([T]) along t and n and multiplying these 
by T and N, respectively, gives ds{aT) = 0, thus T = f = f^dsZ. Since T = 9tX • 5^X, f^d^z = cor^d.cj), and the 
azimuthal slope must also vanish everywhere or nowhere. Observe that the curves never touch the rotational axis 
unless they are coplanar with it, with T = ds4> = C2 = 0. 

Matching sine and cosine terms in the projection of ([T]) along either x or y leads to two more equations, one of 
which has already been satisfied. The remaining equation and constraint equation may now be written in two ways. 
Independent elimination of (f) and z gives the two pairs of equations 
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The limit (ci,9sz) — is problematic for equation pair ([5| and ([6]), and the limit (c2,9s0) —J- is problematic for 
equation pair ([t]) and (|8|. Subsequent manipulations will correspond to multiplications by zero and infinity in these 
limits, so both sets of equations must be retained and examined together. 

Multiply ([5]) by |^ and ([t]) by use ^ and ([8| to rewrite the squares of these terms, then multiply by dsZ or 

rds(t> integrate to get 
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Squaring these and using (|6| and ([s]) again to eliminate dsZ and rds<p gives 
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{dsvf = 1 
{dsrf = 1 
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(9a) 
(9b) 
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The limit (ci, C2) — > is the straight radial line r ~ ±s with (z, (f) constant. Taking C2 
equation [lOj for strings coplanar with the axis of rotation, while taking ci 
[6] for strings in the plane perpendicular to the axis of rotation 



in (10a) leads to Mack's 



in (10b) leads to Fusco's equation 
At first glance, these limits only give sensical 
results if they are taken for the correct equation. However, the products of integration constants in the denominators 
actually represent the same physical quantity, and equations (10) are, in fact, the same equation. Let's thus define 
C3 = C1C311 = C2C3_L, and take stock of our constants of integration. They may be expressed in terms of three vectors, 
namely the velocity, stress vector, and unit axis of rotation: 

(11) 
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These expressions allow recovery of the helical solutions after setting dgV = in ( 10 ). The first two constants represent 
measures of the axial slope and azimuthal slope, while the third is suggestive of an energy. One could divide these 
by the mass density and suitable powers of frequency, and then normalize by a relevant length scale such as the total 
length of string involved in the situation of interest. Presently, there is no such relevant scale, and it seems most 
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sensible to use C3 to provide the length in the normalizing factor. Define the nondimensional arc length s = s ^2^) 

1 

and two parameters a = ci^ and ^ = 02(^1^^ , then substitute 

± dsu = [u^ - 2u^ + (1 - a^)u - /J^] ^ 



u = r ^— in ( 10 ) and fiddle a bit to obtain 



1 

From this, one may write an implicit solution in terms of an elliptic integral 



Expressions for 
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(15) 
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u{S) [1 - u{~s)] 
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If the physical assumptions a > and /i > are made, then relations (11) and (13) imply C3 > and C3 > cf. 
Thus, the cubic coefficients in (14 1 or (15) are such that — 2 < 0, 1 — > 0, and — < 0. Use of "Descartes's 
rule of signs" shows that there are no real negative roots of the cubic — there can be three or one real positive roots, 
depending on the relative magnitudes of the parameters. If /3 = 0, solutions are coplanar with the axis of rotation, 
one root is zero and the other two are 1 ± a. If a = 0, solutions are perpendicular to the axis of rotation, and there 
will be three real roots when < one real root otherwise. In the general case, it is a tedious but straightforward 
process to find all three roots in terms of the slope parameters a and /3, and explicitly reconstruct the parameters in 
terms of the roots. 

For non-helical solutions to exist, both the cubic and u must be non-negative. With three real roots, there are two 
branches of solutions with non-constant radius: an inner branch where u is larger than the smallest root but smaller 
than the other two, and an outer branch where u is larger than all three roots. With one real root, there is only an 
outer branch, where u is larger than the real root. Note, however, that the outer branch of the cubic will correspond 
to a configuration that is concave outward from the axis of rotation, and thus in a state of compressive stress. This 
may be easily seen algebraically for the coplanar case, where the outer branch corresponds tou>l + |a|>l, which 
implies through ( 13 ) that the stress is compressive. A compressive stress in ([T]) corresponds to inherently unstable 
dynamics that arc essentially impossible to observe in real life. So the important physical solutions are those of the 
inner branch, which only exists for sufficiently small relative values of /3 such that the cubic discriminant is non- 
negative. This bound on one of the parameters was previously noted, through a different argument, in Fusco's study 
of periodic, perpendicular-planar strings [6 . The inner branch coplanar solutions are the only spatial curves that 
touch the axis of rotation. Explicit parametrizations of these curves are known, with z expressed as a function of r 

[281 Eg. 

Much of this behavior is more easily seen in phase portraits for u, obtained from equation (14). Four such portraits, 
displaying different choices of slope parameters, are shown in Figure [T] They consist of a closed inner tensile lobe 
and open outer compressive lobe. Increasing either parameter results in smaller lobes. The lobes merge into the line 
u = 1 and a single parabola {dgu)'^ — 4u when both a and (3 go to zero; this is the straight radial line solution. The 
definition of the s coordinate as arc length constrains all trajectory curves to lie inside this parabola. The curves 
intersect the u axis at the roots of the cubic. Figure 1(a) which looks like a fish, shows perpendicular-planar solutions. 

= 0, Z?'^ = 0), the roots are (0, 1, 1) and 
Increasing a as in Figure 1(b) shrinks the lobes, the inner lobe shifting 
towards the dgU axis; the two smallest inner lobe curves have disappeared. Figure 
inner lobe curves touch the dsu axis, and do so with a generically non zero slope ds 

For (a — 1, P'^ = 0), the roots are (0, 0, 2). Increasing /? as in Figure [l(d) shrinks the lobes, the inner lobe shifting 



The inner lobe curves do not touch the dgu axis for nonzero /3 
for a = 0,/3^ = ^, the roots are (|, |, | 



1(c) shows coplanar solutions. The 
although dgU — when u — 0. 



off of the dgU axis; the two smallest inner lobe curves have disappeared. 

A few spatial curves are shown in Figure [2] They resemble jump ropes for small /3 and fiower petals for small a, 
and may be extended indefinitely to any length. Any curve's chirality may be reversed to produce another solution, 
and rotation may occur in either sense, as changing the sign of lu does not affect the equations for r, (j>, and z. 

Consider applying boundary conditions to an end of the string, at a given radius corresponding to a nondimensional 
squared radius uq that lies within one of the ranges allowed by the cubic. Here, and are the values 
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FIG. 1. Phase portraits for u, obtained from equation (141. All curves lie inside the parabola {dju) — Au and intersect the 



+ (1 - a^)u - 



0. 



Increasing either a or /3 shrinks the traj ectories 
(a) The perpendicular-planar case a = 0, = (0.01,0.1,0.3,0.6,0.9). The curves do not touch the dgu axis, (b) a — 0.6 

= (0.01,0.1,0.3,0.6,0.9) ,[(c)]The coplanar case a = (0.01,0.1,0.3,0.6,0.9) , ^/3^ = 0. The curves touch the dsu axis 
a = (0.01,0.1,0.3,0.6,0.9), f /3 = 0.6. 



(d) 



of the axial and azimuthal slopes, respectively, and one must apply a nondimensional stress ^ — 1 — uq along the 
tangent. Requiring periodicity of X restricts the problem to the perpendicular plane [B], while requiring a vanishing 
radius yields only the coplanar solutions. It is likely that a vanishing stress condition, such as would be found at a 
free end, is similarly restrictive for the corresponding problem in the presence of gravity. Imposition of one or more 
of these conditions at the outset of considering this problem, or an augmented form of it, may thus lead to neglect of 
three-dimensional solutions. 

This short study may be expanded in several ways. A logical next step, especially if one is concerned with industrial 
yarn problems, is to add a constant tangential velocity Tq to the string so that 9tX — >■ (?tX-f Tgt and material elements 
move along the rigid shape. From Mack's analysis [lOj . it appears that this will add terms proportional to inside 
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FIG. 2. A menagerie of non-planar jum prop ica. 



arrowed lines are shown to guide the eye. (a) 
a = (0.01,0.1,0.5), f = 0.5 



Vertical dotted lines represent the rotational a xis, and orthogonal 

0.01 



Three views of a 



(0.01,0.1,0.5,0.9), f 



(b) 



Three views of 



the rightmost parentheses in equations (10). This wiU not change the order of equation (14), but should add new 
solutions via a third parameter in the equation that represents the relative strength of Coriolis forces. This string 
problem contrasts with those that occur in inertial, non-rotating frames, where adding tangential motion leads only 
to centripetal forces. Unlike Coriolis forces, centripetal forces are balanced merely by a change in the stress |39j . 
which does not modify the shapes of the solutions, but might stabilize some portion of those that would otherwise 
bear compressive stresses. 

The addition of a uniform gravitational force along the axis of rotation will likely have the nontrivial effect of 
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destroying the periodicity of r, as the solutions of the corresponding hnearized problem are Bessel functions. An 
interesting set of non-rigid solutions to look for, perhaps using the rigid helical solutions as a starting point, are the 
"subharmonic modes of rotation" observed by Caughey, "in which the fixed end of the chain performs an integral 
number of rotations in the time that the free end takes to make one complete revolution" 

Recent numerical work by Aristoff and Stone 0^ suggests that the addition of air drag to a coplanar solution, 
while keeping the boundary condition r = 0, leads to curvature singularities at the radial maxima. Understanding 
the effects of drag on rotating thin structures has direct relevance to the design of wind and water turbines. It is 
curious that, while the patent for the Darrieus turbine |26| cites a rotating rope's shape as that which will minimize 
the bending moments on a slender blade, the helical shape of this turbine's progeny, the Gorlov turbine [40], was 
chosen for other reasons, apparently unrecognized as another choice that would satisfy the same criterion. 

Finally, those forbidden regions of phase space which correspond to negative u might be explored by interpreting 
Lo as complex, and considering a scattering problem for strings in a cylindrical potential. 
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